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SUPER DECOUPLED LOAD FLOW ME TODS 

ABSTRACT 

The complete patent specification involves the invention of 
six super decoupled loadf low methods. The Fast Super Decoupled 
Load flow (PSDL) method and two variants Transformed Fast 
Decoupled Loadflow (TFDLXB and TFDLBX) methods are known ones and 
included are their novel versions (NFSDL, NTFDLXB, NTFDLBX) . 
These are the best versions of many simple variants with almost 
similar performance* These methods specifically involved the use 
of the following invented technics in the prior art method. 

1. Gain matrices of all the six methods are different and they 
can be determined as described in the specification* 
Obviously the gain matrices of pe-subproblem can be defined 
unsymmetrical . 

2. In all the six methods , rotation angles are restricted to 
the maximum value of -36 degrees from nonlinearity 
considerations. 

3. Slack-start for any decoupled loadflow method for efficiency 

4. Modification of real power mismatches at PV-nodes according 
to relations (12) and (15) of the specifications for the 
methods FSDL, TFDLXB, TFDLBX. 

5. - Computation of angle correction and updating for PV-nodes 

along with the voltage magnitude corrections at PQ-nodes for 
NFSDL, NTFDLXB and NTFDLBX methods. In the back-substitution 
part of the solution of this subproblem involving angle 
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corrections at PV-nodes and voltage magnitude corrections at 
PQ-nodes, skip factor elements corresponding to rows and 
columns of PV-nodes in the calculation of voltage magnitude 
corrections. 

The specification deals only with unadjusted super decoupled 

* 

loadflov solution methods and does not describe any applications. 
However the computer algorithms can be appropriately modified 
for adjustments and/or applications such as state estimation, 
contingency analysis etc. 

Also described and claimed is the invention of two compact 
storage schemes for gain matrices of the FSDL method and two 
efficient procedures for node type switching implementations in 
the FSDL method. 



Dated this 7 th day of November 1993 



Signature of the inventor - s. B.Pat el 
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The following specification particularly describes and ascertains 
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the nature of this invention and the manner in which it is to be 
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This invention relates to the super decoupled methods of 
steady-state loadflov analysis of power system, and computer 
algorithms for carrying out these methods. 

A computer algorithm processes raw information to yield 
useful information* A chemical process processes raw material to 
yield useful material. Useful information is well recognized as 
technological product in the modern age of information 
technology. Therefore, the computer algorithm like chemical 
process yields useful product and it is the useful art. 

INTRODUCTION 

Loadflow in power system studies is the most basic frequen- 
tly performed steady state analysis of an electrical power net-- 
work. The loadflows are performed in system planning, operational 
planning, and operation control* Of various methods proposed over 
the past four decades, Stott's Past Decoupled Loadflow (FDL) [1] 
has gained wide acceptance both for off- and on- line applica- 
tions; However, it is known to suffer poor convergence for 
systems having high R/X ratio branches (discussions in [1] & 

■ 

(2)). 

Dy Liacco and Ramarao (discussion of [2]) as well as Deckmann 
et. al. (3] . have developed circuit transformation schemes to 
avoid such difficulties. The two schemes however do not provide 
consistent improvement in convergence [4 }• Also they are not 

♦ 

general and are suitable where the system has a small number of 
troublesome R/X ratio branches. 

A recently introduced super decoupled approach [5,6,7) is 
more general and reliable. Rotation operators applied to the 
complex node injections and the corresponding admittance values 
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that relate the above to the system state variables, transform 
the network equations such that branch admittances appear to be 
almost entirely reactive. Thus, better decoupling is realized. 
However, the super decoupled algorithms are not efficient and are 
not very convenient in contingency analysis. A recent modifica- 
tion to the FDL method for networks with high R/X ratio branches 
[4] is in fact variation of the super decoupled approach with a 
weakness that it takes higher number of Iterations for normal 
cases. 

General-purpose version of the Fast Decoupled Loadflow 
(GFDL) method( 8 } proposes simple modifications to the classical 
FDL method. The convergence is much improved in the presence of 
large R/X ratio lines. The two versions of FDL method are put in 
a better theoretical framework by Monticelli et. al.[9). The GFDL 
method is nothing but an experimental investigation of the 
original observation made about the behaviour of the FDL method 
by this inventor as back as 1985 in an unaccepted research paper 
communicated to IEEE (New York). 

The critically coupled loadflow methods are based on the use 
of the overlap update rule [10]. The disadvantage of the methods 
is that they involve the solution of a larger number of equa~ 

♦ 

tions at each iteration. Moreover the tests conducted for the 
paper [10] avoided the increase of R/X ratio of branches connected 
to the PV-nodes. Therefore, conclusions of the paper are not 
based on exhautive testing. Test show that the invented six 
versions are the best loadflow methods. However, simply modified 
versions including hybrids of the six invented methods can also 
have closely similar performance for any given system. 
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The latest Novel Decoupled Loadflow Method [11] can simply be 
derived from published[12] and lately claimed invetion of this 
patent specification. . In fact the idea involved in the method of 
[11] was revealed by this inventor in extensive communication for 
about four years with Institution of Electrical Engineers of 
England. Also the idea was investigated along with the modifica- 
tion of the real power mismatches at PV-nodes and only the best 
method was reported in reference [12]. Moreover the unsymmetrical 
gain matrix definition for the Pesubproblem is obvious from the 
definition of the factor (K) used to modify real power mismatches 
of the PV-nodes [ 12 ] . However the method of [11] is not reliable 
for systems having PV-nodes and it has no flexibility of tuning 
rotation angle as was originally known to this inventor. Possibly 
this may be the reason that the author of [11] did not comapre 
his method with the Fast Super Decoupled Loadflow [12] and 

■ 

Transformation based Fast Decoupled Loadflow [13] methods. 

Keywords : Power systems, Loadflow, Matrices, solution of 

simultaneous equations. 

The invention will now be described, wherein the following 
symbols are used : 

Y «G +jB ; (p-q)-th element of nodal admittance 
pq pq pq matrix formed excluding shunts 

y = g +jb : total shunt admittance at node p 
p p p 

V =e +jf =»V l& ; complex voltage at node p 
P P P P P 

tfi $ & y 2 voltage angle, magnitude corrections 
P P 
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£e , : real, imeginary component corrections of voltage 

P P 

P + jQ : net nodal injected power calculated 

P P 



♦ * 



+j : nodal power residue (mismatch) 
P P 

RP +jRQ : modified nodal power residue 
P P 

m : number of PQ-nodes 



[ ] : indicates enclosed variables to be vector or matrix 



SUPER DECOUPLED LOADFLOW (SDL) METHOD { The prior art) 



Decoupled Load flow (FDL) method. The one based on the XB-version 
was developed in 1985 [7] and the other based on the BX-version 
is obvious from the BX-version of the FDL method developed in 



admittance (appendix) . 
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k : number of PV-nodes . -^v. ■ 



n=m+k+l : total number of system nodes 

q>p : q is the node adjacent to node p excluding the case q =P Jx£$M 
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There are two versions of the Super Decoupled Loadf low iv.^lS 
method. These are based on the XB- and BX- versions of the Fast " 'V 
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1989 [8]. These methods involve the iterative solution of system 
of equations (1) and (2) . Unlike FDL model the gain matrices in 
its, transformed version , are unsymmetrical because mostly 
different rotations are required to be applied at the terminal 
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nodes of a branch. Symmetrical gain matrices can be obtained by J'4$B 



the. Haley and Ayres technique [7] of applying average of the 



rotations at the terminal nodes of a branch to the branch \$?M'*, 
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SDL based on the XB-version of the PDL (SDLXB) 

The SDLXB-version is described by the following equations 
(1) to (10) . 

[RP] = [Y'J [ ^9] (1) 
[RQ] = [Y") [ fcVJ (2) 

Wherein each element of [RP] and [RQ] at PQ-nodes are given by 
(3) and (4) respectively. Whereas an element of [RP] at PV-nodes 
is given by (5). 

RP * ( bP Cosp + t& Sin* )/V (3) 
P P P P P P 

RQ » (-^P Sinp + kQ CosJJ )/V (4) 
P P P P P P 

RP = &P /V (5) 
P P P 

Trigonometric functions in (3) and (4) and elements of [ Y' ] and 
[Y") are given by (6) , (7) and, <8), (9) and (10). 



P PP ^ PP PP 



2 2 

COBfi = -B / /G "f B (6) 



y"2 2 
G + B (7) 
PP PP 

Y' = -1 / X and Y" » -B (8) 

pq pq pq pq 

= Z -Y' and Y» = -2b'+£-Y" <*> 

pp t*r pq pq p n+r pq 

b' =■ b Cos* or o b (10) 
P P P P 

Where X is the transformed branch reactance defined in appendix 
by eqat¥on (41) and B is the corresponding transformed element 

pq 

of the susceptance matrix. 



•V. • 
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SDL based on the BX-version of the FDL (SDLBX) 

The SDLBX-version differs from the SDLXB-version only in 
relation (8) as given by relation (11). The SDLBX-version 
consists of relations (l) to (7) , (11), (9) and (10). 



INVENTED SUPER DECOUPLED LOADFLOW METHODS 

Six invented super decoupled loadflow models are described 
in this section. Fast Super Decoupled and two versions of .the 
Transformed Fast Decoupled models are Known ones [12,13] and 



Fast Super Decuopled Loadflow (FSDL) 



system of equations (1) and (2). The model cosists of relations 
(1) to (4), (6), (7), (12) to (16), (9) and (10). 

RP ° ^/(KV) (12) 
P P P P 

Elements of [Y'J and [Y N > and the multiplier K in (12) are given 



by 



Y' « -Y for branch r/x ratio < 2.0 

pq pq 



Y n = -y for branch r/x ratio £2.0 

pq pq 

« -(B +0.9(Y -B )) for branch r/x ratio > 2.0 (14) 

pq pq pq 




C£l ZS+if . * - . • ■ •• 



. .... » / 



* ■ 

•■• .v.- 



*' = -B and Y« - -l / X (li) 

pq pq pq pq 
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included are their novel versions. 



**l • » * *. 

The FSDL method involves the iterative solution of the 
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-{B +0.9(Y -B )) for branch r/x ratio > 2.0 

pq pq pq 

--^»q Bp, for branches connected between 

* two FV-nodes or a PV-node and 

the slack-node (13) 



. . - ft 



V 
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i 



K - Absolute (B /Y # ) (15) 
P pp pp 

KKp - Absolute (Y'pp/Y"pp) (16) 



Branch admittance magnitude in (13) and (14) is of the same 

algebraic sign as its susceptance. Elements of the two gain 

matrices differ in that diagonal elements of [Y M ] additionally 

contain the b' values given by equation (10) and in respect of 

elements corresponding to branches connected between two PV-nodes 

or a FV-node and the slack-node. The factor 0.9 in the relations 

(13) and (14) is tuned only for rotation limited to -36 degrees. 

With different Limiting Rotation Angle (LRA) , it needs to be tuned 

again. In two simple variations of the FSOL method, one is to 

make Y" = Y' and the other is to make Y' « Y" 

pq pq pq pq 

■ 

Transformation based Fast Decoupled Load flow (TFDL) 

■ 

The TFDL model is similar to the FSDL model. They differ 
only in the definition of the gain matrices. 
The TFDL ( XB- vers ion ) (TFDLXB) Load flow 

The TFDLXB- vers ion consists of relations (1) to (4), (12), 
(15), (6), (7), (17), (18), (19), (9) and (10). 



-Bpq for branches connected between two PV- 

nodes or a PV-node and the slack-node 



Y'pq - 

-1 / X for all other branches (17) 

L pq 

Y" « -B (18) 

pq . pq 

KK = Absolute (Y f /(2.1/X )) (19) 
P PP f*f pq 

8 



* 
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The TFDL(BX- vers ion) (TFDLBX) Loadflow 

The TFDLBX -vers ion consists of relations (1) to (4), 
(15), (6), (7), (20), (21), (22), (9) and (10). 



(12), 



Y'pq 



pq 



-Bpq 



-B 



-1 / X 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



for all other branches 



pq 



(20) 



(21) 



• • * 
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KK = Absolute (Y' /( 21b )) 

p pp i-»r pq 



(22) 



Where X is the transformed branch reactance defined in appendix 

pq 

by eqation (41) and B is the corresponding transformed element 

pq 

■ 

of the susceptance matrix. 



The factor KK is to be multiplied to the real power mis- 

P 

match at a PV-node switched to PQ-type in node type switching 
implementation. From general considerations, K and KK are 



P 



P 



restricted to the minimum value of 0.75 for FSDL, TFDLXB and 
TFDLBX methods. However they can be tuned for the best possible 
convergence. for any given system. 

• .... 

In all the SDLXB, SDLBX, FSDL, TFDLXB and TFDLBX models 
[Y') and [Y w J are real, sparse, symmetrical and built only from 
network elements. Since they are constant, they need to be 
factorized once only at the start of the solution. Equations (1) 
and (2) are to be solved repeatedly by forward and backward 
substitutions. * 
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[V] and [Y n ] are of the same dimensions (m+k) * (m+k) when 
only a row/column of the slack-node is excluded and both are 
triangularized using the same ordering regardless of the node* 
types. For a row/ column corresponding to a PV-node excluded in 
[Y"], use a large diagonal to mask out the effects of the off* 
diagonal terms. When the node is switched to the PQ-state the 
row/column is reactivated by removing the large diagonal. This 
technique is especially useful in the treatment of PV-nodes in 
the matrix [Y»]. 

It is invented to make this technique efficient while 
solving (2) for by skipping all PV-nodes and factor elements 
with indices corresponding to PV-nodes. In other words efficiency 
can be realized by skipping operations on rows/columns 
corresponding to PV-nodes in the forward-backward solution of (2) 
for (^V. This has been implemented and the time saving of about 
4% of the total solution time (including input/ output) could be 
realized in 14-14 iterations required to solve 116-node system 
with the uniform R-scale factor 4 applied. The time saving has 
been assessed on PC-XT . It should be noted that the same indexing 
and addressing information can be used for the storage of both 
the matrices as they are of the same dimensions and sparsity 
structure . 



Novel Past Super Decoupled Load flow (NFSDL) 



following system of equations (23) and (24). 

[RPJ - [Yl] ( fifi] (23) 



ERQ(PQ-nodes)! « fol [kv(PQ-nodes)l 
RP( PV-node sjj I J \&9 ( FV-nodes)J 



. > ... 



The NFSDL method involves the iterative solution of the 



,". \ . 



.•",V".' 

(24) fr'^Sfcf? 

.».*,% r C*V.\SV 
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Wherein each element of [RP] and [RQ] at PQ-nodes are given by 
(3) and (4) respectively. Whereas an element of (RP) at PV-nodes 
is given by (5). In (24) PV-nodes are assumed to be grouped and 
numbered after all PQ-nodes. 
Elements of (Yl) and (Y2) are given by : 



f - * 



• ■ » - _ _ • 



pq 



-Y 



pq 

(B +0.9(Y -B )) 

pq pq pq 



-Bpq 



Y2pq = 



Gpq 



V 



for branch r/x ratio £2.0 



for branch r/x ratio > 2.0 



for branches connected between 
two PV-nodes or a PV-node and 
the slack-node (25) 



for branches connected to a PQ-node 
and a PV-node 



-Ylpq for branches connecting two PQ-nodes 



Yl for branches connecting two PV-nodes 

pq 



(26) 



■ I 4 



■ ■ * - - * ■ 



■ •-•!•• 



•. .■• y 
... v ■ 

. % . -• 



Yl 



PP 



Y2pp 



£-yi 
*f*p pq 



2b£-Ylpp 



for PQ-nodes 



for PV-nodes 




(27) 



(28) 



-B 

_ PP 

The NFSDL model consists of relations (23), (24), (3) to 
(7), (25) to (28), and (10). Branch admittance magnitude in 
(25) is of the same algebraic sign as its susceptance. 

Hovel Transformed Fast Decoupled Loadf low (NTFDL) 



• ■ ••i ~ 

m 



The NTFDL model is similar to the NFSDL model. They differ 

only in the defination of gain matrices. 

11 
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The NTFDL ( XB-vers ion) (NTFDLXB) Loadflow 

The NTFDLXB model consists of relations (23), (24), (3) to 
(7), (29) to (32) , and (10). 



Ylpq - 



-Bpq 



•1 / X 



pq 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



for all other branches 



(29) 



Y2pq 



Gpq 



B 



pq 



Yl 



PP 



pq 



<H>P pq 



for branches connected to a PQ-node 
and a PV-node 



for branches connecting two PQ-node s 



for branches connecting two PV-nodes (30) 



(31) 



m ♦ 



Y2 



PP 




for PQ-nodes 



for PV-nodes 



(32) 



Where X is the transformed branch reactance defined in appendix 

pq 

by eqation (41) and ST is the corresponding transformed element 

pq 

of the susceptance matrix* 

The NTFDL(BX-version) (NTFDLBX) Loadflow 

The NTFDLBX model consists of relations (23), (24), (3) to 
(7), (33) to (36), and (10). 



Ylpq 



-Bpq 



for branches connected between two PV- 
nodes or a PV-node and the slack-node 



-B 



for all other branches 



pq 



(33) 
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Y2pq 



Yl 



PP 



Y2 



PP 



Gpq 



1/ X 



pq 



pq 



X-yi 
i-*p pq 



for branches connected to a PQ-node 
and a FV-node 

for branches connecting two PQ-nodes 

for branches connecting two PV-nodes (34) 



(35) 



2b' + Z--V X 
P pq 



-B 



PP 



or = 2b'+i/X for PQ-nodes 
P PP 

for PV-nodes (36) 



In the NFSDL NTPDLXB and NTFDLBX methods [Yl) and [Y2] are 
real, sparse, symmetrical and built only from network elements* 
Since they are constant, they need to be factorized once only at 
the start of the solution. Equations (23) and (24) are to be 
solved repeatedly by forward and backward substitutions. 

[Yl] and [Y2] are of the same dimensions (m+k)*(m+k) when 
only a row/column of the slack-node is excluded and both are 
triangularized using the same ordering regardless of the node- 
types. It should be noted that the same indexing and addressing 
information can be used for the storage of both the matrices as 
they are of the same dimensions and sparsity structure. Unlike 
[V], It is to be noted that all the PV-nodes are also active in 
[Y2]. Therefore, the novel methods would pose a difficult problem 
of .node-type switching. 

fi is restricted to the maximum of -36 degrees from 
P 

nonlinearity considerations for FSDL, TFDLXB, TFDLBX, NFSDL, 

NTFDLXB and NTFDLBX methods. However it can be tuned for the 

best possible convergence for any given system and it can also be 

vertually unrstricted -90 degrees. 
■ 13 
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Iteration Scheme 

The basic iterative scheme for solving the proposed models 
is to solve for [Ifi] to update [e] and then solve (2) or (24), 
This is the block Gauss-Seidel approach. The scheme is block- 
successive which imparts increased stability to the solution 

process. This in turn improves convergence and increases 

« 

reliability* Cycling * behaviour in the iterative process of the 
Gtandared iteration scheme of [l] can be avoided by simply using 
different convergence tolerances for the real and reactive power 
mismatches as suggested in [9). The cycling behaviour and its 
remedy as above were originally observed by this inventor as back 
as in 1985 in an unaccepted paper submitted to IEEE and prepared 
from research conducted at the university of Roorkee. However, 
algorithms are given for strictly successive iteration scheme of 

[•]. 

APPENDIX 

Transformation of Branch Admittance 

The branch admittance transforamatlon for symmetrical gain 
matrices of the TFDLXB, TFDLBX, NTFDLXB and NTFDLBX methods is 

* 

given by the following steps : 



2. Compute the average of rotations at the terminal nodes (p 
and q)of a branch 



.. • < v,v. 




• 'v.v'y.Vf.- 
?.~ i • • 

1. Compute : = arctan (G /B ) and 



PP PP PP 

(37) 



fi m arctan (G /B ) 

- ~v ™ . % • • ■ • § - fc 



qq qq qq 



* mm * , » 



ft = (0 + 0 ) / 2 (38) 
av pp qq 



»"■»•• m.tm 

■• .V. 
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3. Compare p with the Limiting Rotation Angle (LRA) and let 

av 



fi to be the smaller of the two : 
av 



p « minimum (0 , LRA) 
av av 



(39) 



V s • 



*.•■"• /• % 
•• . «-.- » 
<■ « 



4. Compute transformed pq-th element of the admittance matrix : 



< # » 
■ ■ ^ » • • ■ 

■ J ' ■ 



G +jB = (Cosp +jsin0 ) (G +jB ) 
pq pq av av pq pq 



(40) 



. *J 



5. Note that the transformed branch reactance is : 



•v. 



2 2 

£ = B / (G +B ) 
pq pq pq pq 



(41) 



and similarly 



% *m * * * 



* •* ■ " * V 



2 2 

I = B / (G +B ) 
PP PP PP PP 



(42) 



• ' * V 

v.f.v 

<-;•.. .*. 



algorithm- l (The prior art) 

* • 

Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (19, IV) Iterative Scheme : 

(a) Read system data and assign an initial approximate solution. 
If better solution estimate is not available, set specified 
voltage magnitudes at PV-nodes and 1.0 p.u. voltage 
magnitudes at PQ-nodes. Set all the node angles equal to 
that of the Slack-node angle. This is referred to as the 
flat-start (11. 

(b) Initialize iteration counts ITRP»ITRQ»r=0 

♦ 

(c) Form nodal admittance matrix 

15 
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(d) Form (m+k) by (m+k) size matrices [Y'] and [Y n ] of (1) and 
(2) respectively each in a compact storage exploiting 
sparsity. 

(i) In case of SDLXB-method, the matrices are formed usig 

* • 

relations (8) , (9) and (10). 

(ii) In case of SDLBX -method, the matrices are formed using 
relations (11) , (9) and (10). 

m 

(e) In [¥"] matrix, replace diagonal elements corresponding to PV 
-nodes by very large value (say 10. 0**10). In case [Y w ] is of 

-*■.■.♦**■-•.-/*:» 

dimension (m by m),this is not regired to be performed. -feMC 
Factorize [ Y' ] and [Y") using the same ordering regardless 
of the node-types. In case (Y"J is of dimension (m by m) , it 
is factorized using different ordering than that of [Y']. 

r r 

(f) Compute residues £P (W" and PV-nodes) and bQ (at PQ- 



(h) Solve (1) for voltage angle corrections and update voltage 
angles using eqn. (43). 

« 

(i) Increament the iteration count ITRP=ITRP+1 and 
r=(ITRP+ITRQ)/2. 

r r 
(j) Compute residues J^P (PQ- and PV-nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( £. ) , 

proceed to step (n) . Otherwise follow the next step. 

16 



. V «A « 



re* ' r • 
••••». •■■ '•• 
* / » ■ * 
.•...\ 

..•'.>.' 



1 ccc* 



nodes only) . If all are less than the tolerance ( £. ) , 



. .".tV 



proceed to step (n) . Otherwise follow the next step. r$&& 



» • * > 

;;■> 



(g) Compute the vector of modified residues [RP] using (3) for 

PQ-nodes and using (5) for PV-nodes. { ^§$0 




■ . * -«'•* a. ■ 

• • * * ■ ' 
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(k) Compute the vector of modified residues [RQ] using (4) for 
PQ-nodes only. 

(1) Solve (2) for voltage magnitude corrections and update 
voltage magnitudes of PQ-nodes only using eqn.(44). 

(m) Increament the iteration count ITRQ=ITRQ+1 and 
r=(ITRP+ITRQ)/2. Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

The SDLXB and SDLBX algorithms differ only in step-d defining 
gain matrices in the above solution steps. Fig.l is the flowchart 
of algorithm-1. 

ALGORITHM- 2 (containing invented steps) 

Fast Super Decoupled Load flow Algorithms (solution Steps) 
Strictly Successive (18, IV) Iterative Scheme | 

(a) Read system data and assign an initial approximate solution. 
If better solution estimate is not available, set all the 
nodes voltage magnitudes and angles equal to those of the 
Slack-node. This is referred to as the slack-start. 

(b) Initialize iteration counts ITRP=»ITRQ=r=0 
. (c) Form nodal admittance matrix 

(d) Form (m+k) by (m+k) size matrices [Y'J and [Y"] of (1) and 
(2) respectively each in a compact storage exploiting 
spars ity. 

(i) In case of FSDL-method, the matrices are formed using 

relations (13), (14), (9) and (10). 

17 
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(ii) In case of TFDLXB-method, the matrices are formed usig 
relations (17), (18), (9) and (10). 

(iii) In case of TFDLBX-method, the matrices are formed 
using relations (20) , (21), (9) and (10). 

(e) In IY W ] matrix, replace diagonal elements corresponding to PV 
-nodes by very large value (say 10. 0**10). In case [Y") is of 
dimension (m by m) , this is not reqired to be done. 
Factorize (Y' ] and (Y"] using the same ordering regardless 
of the node-types. In case [Y M ] is of dimension (m by m) , it 
is f actor ized using different ordering than that of (Y'J. 



r r 
(f) Compute residues (*Q- and PV-nodes) and (at PQ- 



(g) Compute the vector of modified residues [RP] using (3) for 



angle) from nonlinearity considerations. It can be even 
unrestricted for any given system (-90 degree). 



P P P 

(i) Set the voltage magnitudes of PV-nodes equal to the 
specified values. This is required to be done only in the 
first iteration when slack-start is used. Increament the 

m 

iteration count ITRP=»ITRP+1 and r=(ri?*P4- ZTR$)/2-> 

IV 



•* ** *• 

• ♦ * »*- " * ■ t 

.«. '.v.a. 



•V-V *. • 

r » X - ■ _ • • 

• . — " ?>\. . 
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'.".V » .v. /.'Ay 
»;•■•.-. V.V.W. 

"* . -V V/ , ■ .•>• 
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nodes only). If all are less than the tolerance ( ^), 
proceed to step (n) . Otherwise follow the next step. 



PQ-nodes and using (12) and (15) for PV-nodes. In (3) £ is 

P :■<•'■.'« \J*;V 

restricted to the maximum, of -36 degrees (or any other 



- • ■ . J 



v.: p • - 

(h) Solve (1) for voltage angle corrections and update voltage 
angles using relation (43). 



v.- 

.»•.. ■ . • . 

.-V, .-. -v . 



r (r-X) r 

e = e +1^ (43) 



•«! •:*•* . 
. ■ .■. . 
..... 

• * * .* 
»* ■ , ■ * ■ 



*•"•*•.{».*• * 

a..** r - 



- . - 




2107388 



r r 
(j) Compute residues f^P (PQ- and PV-nodes) and (at PQ- 

nodes only). If all are less than the tolerance ( £ ) , 

proceed to step (n) • Otherwise follow the next step. 

(k) Compute the vector of modified residues [RQ] using (4). In 

(4) 0 is restricted to the maximum of -36 degrees (or .any 
P 

other angle) from nonlinearity considerations. It can be 
even unrestricted for any given system (-90 degree)* 

(1) Solve (2) for voltage magnitude corrections and update 
voltage magnitudes of PQ-nodes only as given by relation 
(44). While solving (2), skip all the rows and columns 
corresponding to PV-nodes to increase the efficiency. 

r <r-l) r 

V = V + bV (44) 

P P P 

(m) Increament the iteration count ITRQ=»ITRQ+1 and 
r= ( ITRP+ITRQ)/2 . Proceed to step (f) 

■ 

(n) Calculate line flows and output the desired results. 

All the above steps can also be the parts of the classical 
(18, IV) iteration scheme of reference [1J. This algorithm is for 
solving loadflow problem formulated in polar coordinates. The 
solution of the loadflow model formulated in rectangular 
coordinates also involve similar steps with minor modifications. 
Gain matrices of the invented methods are the same for both 
polar- and rectangular-coordinate formulations of the loadflow 
problem. Pig .2 is the flowchart of algorithm-2. The FSDL-method, 
TFDLXB-method and TFDLBX-method differ only in step-d defining 
gain matrices in the above algorithm-2. 
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ALGORITHM- 3 (invclude invented steps) 

Novel Fast Super Decoupled Loadflow Algorithms (solution Steps) 
Strictly Successive (18, IV) Iterative Scheme | 

(a) Assign an initial approximate solution. If better solution 
estimate is not available, set all the ' nodes voltage 
magnitudes and angles equal to those of the Slack-node. This 
is referred to as the slack-start. 

(e) Initialize iteration counts ITRP=ITRQ=r=0 
(c) Form nodal admittance matrix 



(d) 



Form (m+k) by (m+k) size matrices [Y1J and [Y2] of (23) and 
(24) respectively each in a compact storage exploiting 
spars ity. 

(iv) In case of NFSDL-method, the matrices are formed using 
relations (25) to (28), and (10). 



(v) In case of NTFDLXB-me thod , the matrices 
using relations (29) to (32) and (10) . 

(vi) In case of NTFDLBX-method , the matrices 
using relations (33) to (36) and (10). 



are formed 



are formed 



(e) Factor! ze [Yl] and [Y2] using the same ordering regardless 
* of the node-types. 

r r 

(f) Compute residues &P (PQ- and PV-nodes) and (at PQ- 

nodes only). If all are less than the tolerance ( £_), 
proceed to step (n) . Otherwise follow the next step. 
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Compute the vector of modified residues (RP) using (3) for 

PQ-nodes and using (5) for PV-nodes. In (3) is restricted 

P 

to the maximum of -36 degrees (or any other angle) from 
nonlinearity considerations* It can be even unrestricted for 
any given system (-90 degrees). 

Solve (23) for voltage angle corrections and update voltage 
angles using relation (43). 

Set the voltage magnitudes of PV-nodes equal to the 
specified values. This is required to be done only in the 
first iteration when slack-start is used. Increament the 
iteration count ITRP»ITRP+1 and r~(ITRP+lTRQ)/2. 

r r 
Compute residues (PQ- and PV-nodes) and (at PQ- 

nodes only) . If all are less than the tolerance ( L ) , 
proceed to step (n) • Otherwise follow the next step. 

Compute the vector of modified residues [RQJ using (4) for 

PQ-nodes. Compute [RP] for PV-nodes using (5). In (4) p is 

P 

restricted to the maximum value of -36* degrees (or any other 
angle) from nonlinearity considerations. It can be even 
unrestricted for any given system (-90 degrees) • 

Solve (24) for voltage magnitude corrections at PQ-nodes 
and voltage angle corrections at PV-nodes. Update voltage 
magnitudes of PQ-nodes using relation (44) and angles of PV- 
nodes using relation (43) • In the back substitution part of 
the solution of (24), skip all the PV-nodes in the 

calculation of voltage-magnitude corrections of PQ-nodes. 

• » 

21 
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(m) Increament the iteration count ITRQ-ITR0+1 and 
r=»(ITRP+iTRQ)/2. Proceed to step (f) 

(n) Calculate line flows and output the desired results. 

All the above steps can also be the parts of the classical* 
(ie ( IV) iteration scheme of reference [1]. This algorithm is for 

« 

solving loadflov problem formulated in polar coordinates* The 
solution of the loadflow model formulated in rectangular 
coordinates also involve similar steps with minor modifications* 
Gain matrices of the invented methods are the same for both 
polar* and rectangular-coordinate formulations of the loadflow 
problem. The NFSDL-method, NTFDLXB-method and NTFDLBX -method 
differ only in step-d defining gain matrices in the above 
algorithm-3. Fig. 3 is the flow-chart of algorithm-3. 

ALGORITHMS using GLOBAL CORRECTIONS 

■ 

The algorithms-l # -2 and -3 of above involve increamental (or 
local) corrections. All the above algorithms can be organised to 
produce corrections to the initial estimate solution. It involves 
storage of the vectors of modified residues and replacing the 
relations (3) , (4) , (5) , (12) , (43) and (44) respectivly by 
(45) , (46) , (47) , (48) ,(49) and (50). Superscript '0' in relations 
(49) and (50) indicates the initial solution estimate, 
r r r r. (r-1) 

RP - ( COS0 + fcQ Sinfi )/V + RP (45) 
P P P P P P P 

r r r r (r-1) 

RQ « (-^P Sinp + tfl Cosfi )/V RQ (46) 

P P P P P P P 
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r r r (r-l) 

RP = J^P /V + RP 
P P P P 



(47) 



RP 



r r (r-l) 

&P /(K *V ) + RP 

P P P P 



(48) 



• * 



e 



0 r 

e + 

p p 

o r 

p p 



(49) 



(50) 



RECTANGULAR COORDINATE FORMULATIONS OF THE , INVENTED LOADFLOW 
METHODS 

This involves following changes in the equations describing 
the loadflov models formulated in polar coordinates • 

(i) Replace 0 and fcfi respectively by f and &f in equations 

(1) , (23), (24), (43) and (49). 

i 

(ii) Replace V and respectively by e and ^e In equations 

(2) , (24), (44) and (50). 

• * 

(iii) Replace V by e or e in equations (3), (4), (5), {12), 

p p s 

(45), (46), (47) and (48). The subscript 's' indicates the 
Slack-node variable. 

* 

* 

(iv) After calculation of corrections to the imaginary, part of 
complex voltage ( ^f ) of PV-nodes and updating the 
imaginary component (f) of PV-nodes, calculate real 
component by : 



V 



- f 



2 



p( specified) p 



(51) 
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COMPACT STORAGE SCHEMES FOR GAIN MATRICES OF THE FSDL METHOD 

Two compact matrix storage schemes for gain matrices of the 
FSDL method are described in the following. 



1. Identify PQ-nodes in the composite path of all PV-nodes. 



2. Reload and ref actor all row/columns of these PQ-nodes using 



3. When solving for £V using factors of [Y'J : 

(i) all PV-nodes and factor elements with indices 
corresponding to PV-nodes are skipped. Zeroing of reactive 
power mismatches at PV-nodes is not required. 

(ii) for PQ-nodes of the path, separately stored factors as 
in step-2 above are used along with the other factors of 
[V] in forward-backward operations. 



The above procedure is implemented. For 118 -node system, 
there are only 18 PQ-nodes in the composite path of all 53 PV- 
nodes. Bach of the gain matrices [Y'] and [Y«J is having 117- 
row/columns for 118-node system, and the total of 2 34 -row/columns 
are to be factor! sed and stored for conventional 2 -matrix 
solution. Threfore the proposed second scheme can be said to be 
saving 42% (99/234=0.42) of the factorization effort and core 
memory. This saving rises to the maximum 50% for systems without 
PV-nodes. 

24 
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(1) When network shunts are ignored (FK=0.0), [Y w J becomes the 
submatrix of [Y'J. In this scheme, [Y'J is factorized using 

■ 

optimal ordering regardless of node-types along with 1 the 
following procedure : 
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(2) When network shunts are ignored (FK=0.0), [Y°] becomes the 
submatrix of [¥'). In the single matrix version of the PSDL 
method only [V] is factorized and stored. This is achieved 
exactly by placing PV-nodes after all PQ-nodes and by allowing 
optimal ordering only within each group. While solving (2) for 
£V, all PV-nodes and factor elements with indices corresponding 
to PV-nodes are, skipped* Zeroing of reactive power mismatch at 
PV-nodes is not required. Apperantly, 50% reduction in the matrix 
storage is achieved in the absence of PV-nodes. 

# 

PV-node Q-limit Adjustment for the PSDL method 

When the power flow solution is moderately converged, any Q- 
limit violations should be corrected. Any violated PV-node is 
converted to the PQ-type with the rective generation set at the 
limiting value. The voltage of the converted node is subsequently 
compared to the scheduled value and the node is reconverted to 
the PV-type, if any of the back-off conditions are satisfied. 

The conventional node-type switching approach is more 
suitable for the handling of the PV-node Q- limiting. The 
switching of node-types involves the insertion/delation of 
eqalities into/from (2). This can be achieved efficiently by 
partial matrix refactorization and factor * updating methods. 
Partial Refactorization Method-1 (PR1)[13] can be used to update 
(If") in order to account for changes in the node status (regu- 
lated or nonregulated) . The efficiency of PR1 is inversely 
proportional to the number of node status changes. Two efficient 
techniques of the node-type switching are described in the 
following. 

25 
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(1) When network shunts are ignored (FK=0.0) or accounted by 
alternative ways, a reactivated row/column by the removal of the 
large diagonal corresponding to a FV-node in [Y n ] is identical 
to that in [Y'J. It is this feature of the FSDL method exploited 
to enhance the efficiency of the node type switching using the 
following procedure. This is convenient with the simple variant 
of the FSDL method mentioned before. 

1* find the composite path of PV-nodes switched to PQ-type 

* 

2. find the composite path of the rest of the PV-nodes 

3. ref actor only rows/columns of [Y H J that constitute common 
path of steps 1 and 2. The factored rows/ columns of [ Y' ] in 
the uncommon path of step 1 can be used for the partial 
ref actor izat ion of [Y n ] 

4* while solving (2) , use factors from [Y' J in the uncommon 
part of the path of step 1 along with other factos of [Y"J. 

When all the PV-nodes of the system are switched to PQ-type, 
[Y'J is used for solving (2) and no partial refactorization in 
[Y ,a ] is required. Also when common path to be determined in step- 
3 does not exist, partial refactorization is not required. 

% 

(2) The compact storage scheme (1) for the FSDL method described 
above inherently provides for efficient implimentation of node- 
type switching used for FV-node Q-limit adjustment. The procedure 
constitutes the following steps. 

1. PQ-nodes in the composite path of all PV-nodes are known and 
their rows/columns are factored and stored separately from 
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2. Identify PQ-nodes in the composite path of current PV-nodes. 

3. Identify PQ-nodes which are not common in steps 1 and 2. 

4. Carry out Reverse Correction Partial Matrix Re factorization 
along the composite path of PQ-nodes of step-3 on separately 
stored factors of step-1. Use large diagonals in 
rows/columns of PQ-nodes of step-3. 

If no PQ-node is identified in step-3, Reverse correction Partial 
matrix refactorization of step-4 is not required to be performed. 
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I CLAIM : 

(a) To have invented three Super Decoupled Loadflow methods 
(FSDL, TFDLXB and TFDLBX) both in Polar coordinate and 
Rectangular coordinate 3 formulations, their Hovel versions 
(NPSDL, NTFDLXB and NTFDLBX) , and their all possible hybrid 
cobinations. Also I claim to have invented the technique of 
producing Global corrections as given by relations (45) to 
(50) . 

,b) specifically the invention involves the use of the following 
seven (7) items in the methods of (a): 

■ 

(i) the definitions of the gain matrices of the six-methods 
of (a) corresponding to steps-d in algorithm-2 and 
algorithm-3 , for both polar and rectangular formulations of 
the loadflow problem. Though explicit definitions, of 
asymmetrical gain matrices of Pe-subproblem are not given, 
they are obvious for all the six methods and they are also 
claimed by this inventor as simple modifications. 

(ii) restriction of rotation angle to the maximum of -36 

degrees (any other angle) in relations (3), (4), (45) and (46) 
from nonlinearity considerations, for both polar and 
rectangular coordinate formulations of the loadflow problem. 

(iii) modification of real power mismatch at PV-nodes accor- 
ding to relations (12) and (15) in FSDL, TFDLXB and TFDLBX- 
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• ■ 



methods formulated in both polar and rctangulr coordinates. 
This in general is the technique to modify the known vector 
to ensure the coefficient/gain matrix symmetrical. It can 
be used in all possible computer algorithms (including all 
other approaches of the formulation of the loadf low problem) 
and real-time/on-line/off-line applications. 



solving 
methods. 

(24) , I 



equation 



(iv) computation of angle corrections and updating for PV- 
nodes along with the Voltage magnitude corrections while 

(24) for NFSDL, NTFDLXB and NTFDLBX 
In the back-substitution part of the solution of 
claim to have invented the technique to skip all 
the rows and columns corresponding to PV-nodes for calcula- 
tion of voltage magnitude corrections of PQ-nodes. This also 
applies when the novel methods are formulated in rectangular 
coordinates* 

a 

♦ 

(v) the Slack-start procedure for all the decoupled loadflow 
methods and in particular for the six methods of (a) above. 

(vi) similar decoupled methods as of (a) ban also be used 
for solving simultaneous equations appearing in other 
areas of analysis, operation and control. 

(vii) The six invented methods of unadjusted loadflow 
solution can be used for adjusted solution by adding some 
more steps in the algorithms, state estimation, contingency 
analysis and in variety of advanced power network analysis 
and control. 
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(c) Also I claim to have invented two compact storage schemes 
for gain matrices of the FSDL method and two efficient 
procedures of the node-type switching implementations for 
the FSDL method as described in the above. 



Dated this 7 th day of november 1993. 



V- 



Signature of the Inventor S.B.Patel 
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Fig. t Flow-chart of unadjusted load flow solution aJgorfthm-1 . 
It is the same for SOLXB and SDLBX methods, (flow-chart of the prior art) 
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